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ABSTRACT 

Aims. We present a detailed and uniform study of oxygen abundance from two different oxygen lines at 6158A and 6300A in a large 
sample of solar-type stars. The results are used to check the behaviour of these spectral lines as oxygen abundance indicators and to 
study the evolution of oxygen in thick and thin disk populations of the Galaxy. 

Methods. Equivalent width measurements were carried out for the [Ol] 6158A and Ol 6300A lines. LTE abundances were obtained 
from these two lines in 610 and 535 stars, respectively. We were able to measure oxygen abundance from both indicators in 447 stars, 
enabling us, for the first time, to compare them in a uniform way. Careful error analysis has been performed. 

Results. We found that oxygen abundances derived from the 6158A and 6300A lines agree to within 0. ldex in 58% of the stars in 
our sample, and this result improves for higher signal-to-noise values. We confirm an oxygen enhancement in stars of the thick disk, 
as has also been seen for other ff-elements. The new oxygen abundances confirm previous findings for a progressive linear rise in the 
oxygen-to-iron ratio with a slope equal to 0.78 from solar metallicity to [Fe/H]—1. However, the slope we measured is steeper than 
the one found in previous studies based on the oxygen triplet. Below [Fe/H]=-0.6 our stars show [O/Fe] ratios as high as ~0.8, which 
can be interpreted as evidence for oxygen overproduction in the Galactic thick disk. These high oxygen abundances do not pose a 
problem to chemodynamical models since there is a range of parameters that can accommodate our results. 

Key words, stars: abundances - stars: atmospheres - stars: solar-type - Galaxy: abundances 


1. Introduction 

After hydrogen and helium, oxygen is the most abundant ele¬ 
ment in the Universe. As a matter of fact, oxygen is the only 
chemical element that has a unique production site. This has 
been used by many authors to argue that oxygen is the best tracer 
of galactic chemical evolution. It is produced by alpha particles 
during hydrostatic nucleosynthesis in massive stars. The inter¬ 
stellar medium is enriched with o xygen when massive stars ex¬ 
plode as Type II supernova e, (e.g. IWheeler et al.lll989t iMaedeil 
20091 IStasihska et alJl2012h . 

Despite the importance of deriving precise oxygen abun¬ 
dances in solar-type stars, this is usually not a trivial process. 
The number of spectral lines of atomic oxygen suitable for 
abundance analysis in visual spectra is small. There is a large 
number of molecular lines of OH in the near-UV and near- 
IR, strong enough to be measured in main-sequence metal- 
poor stars down to [Fe /H]—3.5 dex (llsraelian et alJI 19981 1200 lb 
Boesgaar d et al.ll999l) . However, it is extremely difficult to mea¬ 
sure abundances from molecular features as they are affected 
by several physical processes that are not well understood in 
solar-type stars (UV opacities, etc.). The situation with atomic 
lines is no more promising. Serious disagreements exist in both 
solar-type main sequence and evolved sub-giant/giant stars in 
any metallicity range. The well-known spectral lines of atomic 
oxygen suitable for abundance analysis are [Ol] 6300A and the 


Ol triplet at 7771-4A. The weak spectral lines of Ol 6158A 
and [Ol] 6363 A are hardly used owing to their small equiva¬ 
lent widths. The strong lines of the oxygen triplet are strongly 
affected by deviations from local thermodynamic equilibrium 
(LTE) (e.g. ICaffau et al.ll2008l) . which depend on the cross sec¬ 
tion of collisions between the atoms of oxygen and neutral hy¬ 
drogen. The physical description of these collisions is still a 
matter of scientific debate. Moreover, the forbidden line is of¬ 
ten blended with telluric features and also contains a Ni I b lend 
at 6300.40A (e.g. lLamberdll978 AllendePrietoe t al .11200 it) and 
a much weaker CN blend (e.g. Teske et all 2013 ). The Ni line, 
with an uncertain oscillator strength, may contribute up to 60% 
to the equivalent width (EW) of the 6300A feature (see sec¬ 
tion [4]). In addition, a disagreement between the Ol triplet and 
forbidden line abu ndances in metal-poor giants can reach up to 
1 dex (llsraelian et al .1120041) . 

There are many articles dedicated to oxygen measure¬ 
ments in solar-type stars. Several workshops and conferences 
have been held to review the status of oxygen measure¬ 
ments and address unresolved problems related to inconsis¬ 
tent abundances derived from different oxygen lines. Most 
of the previous work on oxygen with large samples of stars 


(e.g. Bensbv et al][2 004: Takeda & Honda 2 005b Ecuvillon et afl 


120061 : Delgado Mena et al. 2010b Petigura & Marcvl 12011 

Ram irez et all2013t) were based on two lines: Ol triplet -7774 A 
and [Ol] 6300A. lEcuvillon et al.l ([2006) have compared near-UV 


Article number, page 1 ofQ3] 














































A&A proofs: manuscript no. oxygen 



[Fe/H] Effective temperature Logfg) 

Fig. 1 . Distribution of the fundamental stellar parameters in the HARPS sample for stars with effective temperatures higher than 5200 K. 


OH lines with the OI triplet and [OI] and found good agreement 
between the [O/H] ratios from forbidden and OH lines, while the 
NLTE (non-local thermodynamic equilibrium) trip l et sho wed a 
systematically lower abundance. iTakeda & Honda! ( 2005 ) have 
used the 6158A OI line to derive oxygen abundances in 160 FGK 
dwarfs and sub-giants. Nevertheless, the quality of their data was 
not high enough to reach conclusions about the reliability of this 
line. 

To all the above-mentioned we must add that the solar oxy¬ 
gen abun dance itself is still uncertain. Since the first determina¬ 
tions o f lLambertl (1 1 978h (log e(0)=8.92) and lAnders & Grevessel 
(Il989h (log e(0)=8.93), the value of solar oxygen has undergone 
several revisions that_yielded progressively lower values. In re¬ 
cent years [Asplund et al.l (f20041) have recommended a value of 
log e(0)=8.66+0.05 based on measurements of [OI], OI, and IR 
OH lines, which was supported by the lSocas-Navarro & Norton! 
( 2007 ) determination (log e(0)=8.63) using the infrared triplet. 
This value was reviewed in 2009 to log e(0)=8.69±0.05 us¬ 
ing new 3D hydrod ynamical models of t he sola r atmosphere 
dAsnlund et al.fl20()9fi . On the other hand lAvresi (i2008t) deter¬ 
mined a significantly larger value (log e(0)=8.81) based on a 
single snapshot of a 3D model for the [OI]6300 line and treat¬ 
ing the log (gf) of Ni blended line as a free parameter. Recently, 
ICaffau et al.l ( 2008 ) presented a new determination of the solar 
photospheric oxygen abundance by analysing different spectral 
atlases of the solar flux and disk-centre intensity, making use of 
the latest generation of C05BOLD 3D solar model atmospheres. 
They have studied the photospheric oxygen abundance by con¬ 
sidering only lines from atomic transitions. Ignoring the role of 
collisions with hydrogen atoms on the NLTE level populations of 
oxygen, they proposed log e(0)=8.76+0.07. They have stressed 
that the measurement of equivalent widths with high precision is 
still an important and open issue. Strangely enough, even for the 
solar high-resolution atlas dKurucz et al.lll984l) . they concluded 
that the placement of the continuum is not trivial, and that the 
influence of many blends gives rise to significant uncertainties. 

Our work presents a complete and uniform study of the oxy¬ 
gen abundances in a large sample of stars from the HARPS sur¬ 
vey, using two different abundance indicators: the high excitation 
line at 6158A and the oxygen forbidden line at 6300A . This is 
the first systematic comparison of these lines in a large sample 
of main-sequence solar-type stars. The paper is organized as fol¬ 
lows: in Sect. 2, we introduce the sample used in this work. The 
method of the chemical abundance determination, analysis and 
uncertainties calculations will be explained in Sect. 3. Sect. 4 is 
dedicated to the comparison between 6158A and 6300A as oxy¬ 
gen abundance indicators, where we provide details about their 
reliability in different ranges of stellar parameters. A discussion 


of the [O/Fe] trends, separation of the different populations of 
stars and comparison with previous works can be found in Sect. 
5. Finally, models of oxygen production in the Galaxy and com¬ 
parison with observations is presented in section 6. 


2. Sample description and stellar parameters 


In the past decade a huge effort has been made to discover ex¬ 
trasolar planets. Several radial velocity surveys have carried out 
a thorough follow-up of solar-type stars over several years, re¬ 
sulting in homogeneous homogeneous sets of high resolution 
spectra and large catalogues of uniform stellar parameters (e.g. 
SWEET-Cat - ISantos etal, 1120131) . Besides the well-known suc¬ 
cess of such projects for discovering extrasolar planets, these ob¬ 
servations give us, for the first time, the opportunity to study the 
atmospheres and the chemical compositions of solar-type stars 
with an unprecedented high-quality, homogeneous and extended 
dataset. 

The sample of stars used in this work is a combination of 
three subs amples from the H ARPS plane t search programme : 
HARPS-1 (iMavor et al.ll2003l) HAR PS-2 dLo Curto et al.ll20Toh 
and HARPS-4 (ISantos et alJ 1201 ll) . These stars are mostly 
slowly rotating chromospherically inactive, unevolved solar-type 
dwarfs and some sub-giants, with spectral types between F2 
and M0. The individual spectra were reduced with the HARPS 
pipeline and then combined using IRAlQ after correcting for its 
radial velocity. The combined spectra for each star have a resolv¬ 
ing power R ~ 115000 and a signal-to-noise (S/N) between ~40 
and ~2000. About 90% of the stars have an S/N higher than 100. 
This high-quality data allowed us to measure the 6158 A line in 
a large sample of stars and investigate its reliability as an oxygen 
abundance indicator. 

We selected 762 stars from the HARPS samples with tem¬ 
peratures higher than 5200K, for which our stellar parameters 
are most precise (Tsantaki et al. 2 0131). S tellar parameters were 
determined by Sou sa et alJ (12008 1201 lalfbh using the same spec¬ 


tra that we studied. The typical precision uncertainties are of 
about 30 K for T e ff, 0.06 dex for log(g), 0.08 km s _1 for and 
0.03 dex for [Fe/H]. Figure Q] shows the distribution of stellar 
parameters for our final sample. 


3. Abundance analysis 

Oxygen abundances were derived from OI6158A and 
[OI]6300A spectral lines. There are two main issues that 

1 IRAF is distributed by National Optical Astronomy Observatories, 
operated by the Association of Universities for Research in Astronomy, 
Inc., under contract with the National Science Foundation, USA. 
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Fig. 2. Spectral regions containing oxygen lines in HD213240, a star from our database with an average S/N = 300 representative for our catalogue. 
HD98281, for which the high S/N allows us to measure lines as small as 1 mA, and for the Kurucz Sun Atlas ( [Kurucz et alJI 19841) . Stellar 
parameters, equivalent width of the oxygen lines (Ni blend contribution removed) and LTE abundances are provided in the plot. We note that the 
provided signal-to-noise ratio is the nominal S/N. 



Fig. 3. Distribution of measured equivalent widths for each spectral line. 


make the automatic EW measurement of these lines prob¬ 
lematic. First, the two lines are placed in complex spectral 
regions close to some strong Fel lines and, in the case of the 
6300A transition, also with telluric lines in its vicinity (Fig. |2}. 
Moreover, the strength of these lines in F- and G-type stars is 
generally smaller than 8 mA (Fig.|3j. Altogether, this makes the 
continuum placement and EW measurement a serious challenge 
for any automatic code. Therefore, EW measurements were 
carried out for both lines in a detailed individual analysis using 
the splot package from IRAF 1 . 

In combined spectra the same atmospheric lines can appear 
at different wavelengths, slightly shifted between them and with 
different intensities. To account for possible blends with telluric 
lines, we made a thorough observation of the spectra in the re¬ 
gion of the forbidden line. This process was done by identifying 
the atmospheric pattern, which is repeated along the region of 
[OI]6300A, in a nearby window free from stellar lines. This pat- 


Article number, page 3 ofED 














































A&A proofs: manuscript no. oxygen 



Fig. 4. [O/H] vs. [Fe/H] for different indicators. The yellow dot corresponds to the solar abundances. 


Table 1. Atomic parameters of the spectral lines together with our mea¬ 
surements of equivalent width and abundance in the Sun. This solar 
value will be used as reference in the present analysis. 



A 

(A) 

log (gf) 

Xlo 

(eV) 

Ref. 

EW 0 

(mA) 

loge(X) 0 

OI 

6158.171 

-0.296 

10.74 

1 

3.7 

8.71 

[OI] 

6300.304 

-9.717 

0.00 

1 

3.4 

8.65 

Nil 

6300.336 

-2.110 

4.27 

2 

2.0 

6.25 


References. (l) ICaffau et all 1200 8ll . (2l l.lohansson et all d2003t) 


tern can be later extrapolated to check if it is superimposed to 
the 6300A line. On the other hand, de-blending was required for 
the 6158A feature in order to account for the presence of the Fel 
line at 6157.5A. The red wing of this Fe line frequently over¬ 
laps with the oxygen line. The spectral windows containing both 
oxygen spectral lines are shown in Figure U For some still un¬ 
known reasons, the nearby line of oxygen at 6156.75A provides 
very different abundances from those based on the 6158A and 
6300A lines. In fact, being even stronger than the line at 6158A, 
it has been avoided in previous studies (e.g. ICaffau et al.l f2008: 
iTakeda & Hondafl20()5l) . _ _ 

The Solar Atlas of Kurucz ( Kurucz et al. 1984) has been 
used to select continuum regions free from spectral lines. 
We found that the continuum windows are 6156.94-6157.10A 
and 6158.80-6158.90A for OI 6158A, and 6298.20-6298.80A, 
6297.15-6297.45A and 6303.00-6303.30A for the [OI] 6300A 
line. We have also selected two stars in our sample, with temper¬ 
atures far from that of the Sun, in order to verify the reliability 
of those regions at different T e ff- This procedure guarantees a 
uniform continuum placement in our targets, whose spectra are 
similar to those of the reference stars. 

Oxygen abundances from two indicators were determined 
according to a standard LTE analysis with the revised version 
of the MOOG2013 spectral synthesis code dSnedenl 19731) . using 
the abfind driver and a grid of Kurucz ATLAS9 plane-parallel 
model atmospheres dKurucz 1993 ). Parameters of spectral lines 
considered in our analysis are listed in TableQ] NLTE corrections 
for OI6158A are negligible (e.g. ICaffau et alj2008l) and were not 
taken into account. On the other hand [OI]6300A is not affected 
by deviations from LTE. The two oxygen lines have very dif¬ 
ferent excitation energies, but very similar strengths in the Solar 


Atlas dKurucz et alll984i) . Once we move towards hotter/cooler 
stars along the main sequence, one of these lines will become 
stronger and the other weaker. Thanks to the wide range of T e ff 
and [Fe/H] of the stars in our sample (Fig. [TJ- we have an inter¬ 
esting opportunity to compare oxygen abundances derived from 
two lines with very different formation physics. We note that, 
for still unknown reasons, these lines provide different oxygen 
abundances in the Sun dCaffau et al.ll2008[) . One of the important 
questions that will be addressed in this article is whether or not 
this difference is seen in other sun-like stars. In general, we want 
to understand which of these lines is more reliable as an oxygen 
abundance indicator and why. 

It is well known that the forbidden line at 6300.30 A is 
severely blended with the Nil 6300.34A line. The contribution 
of the Nil line is negligible at [Fe/H] < -0.8, but becomes very 
important at solar metallicities (see section 0. The EW of the 
Nil line in our stars was estimated using the ew find d river of 
MOOG and Ni abundances from lAdibekvan et alJ (!2012h . calcu¬ 
lated from the same spectra and stellar parameters used in this 
work. The derived EW of the Ni line was subtracted from the 
total EW of the 6300A feature to obtain the contribution of oxy¬ 
gen. 

Many authors use the so-called solar gf-values to carry out 
differential abundance analysis, meaning that they force the os¬ 
cillators strengths of different spectral lines of the same chemi¬ 
cal element to provide a unique abundance value. Here we take 
a different approach. We do not modify the gf-values of the oxy¬ 
gen lines to obtain the same abundance. In this approach we ob¬ 
tain 0.06 dex difference between the solar oxygen derived from 
the 6158A and 6300A lines. We want to investigate whether this 
difference exists in other sun-like stars with a quality of spectra 
similar to the Solar Atlas. The solar abundances listed in TableQ] 
are used as a reference in the analysis. 

Abundances from OI6158A were obtained for 535 stars, 
while 610 stars were analysed using [OI]6300A from the initial 
sample of 762 stars. The final results from both lines are shown 
in Figure [4] These are presented relative to the solar reference 
abundance (TableQ}, derived following the same procedure as 
for the rest of the sample. Tables |2]and[3]provide the abundances 
from the OI6158A and [OI]6300A lines, respectively, together 
with the stellar parameters, the equivalent width of the lines and 
its total uncertainty. 

In Fig. [5] we plot oxygen abundances against stellar param¬ 
eters. Such plots reveal systematic effects that could influence 
our measurements. We note a slight dependence of [O 6300 /H] on 
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Fig. 5. Systematic trends of [O/H] for different indicators with temperature and surface gravity. The solid line is the cubic best fit. 


Table 2. Oxygen abundances from the OI 6158A line 


Star 

T eff 

(K) 

log(g) 
(cm s~ 2 ) 

ft 

(km s' 2 ) 

[Fe/H] 

EW 0/ 

(mA) 

c5EW 0/ 

(mA) 

[0/H] 0/ 

Gal.Pop. 

Sun 

5777 

4.44 

1.00 

0.00 

3.7 

0.3 

0.00 

thin 

HD70889 

6051 

4.49 

1.13 

0.11 

6.9 

0.5 

0.10 

thin 

HD21161 

5923 

4.24 

1.14 

0.09 

7.5 

2.8 

0.18 

thin 

HD37226 

6178 

4.16 

1.61 

-0.12 

7.4 

0.8 

-0.01 

thin 

HD125881 

6036 

4.49 

1.10 

0.06 

5.8 

0.6 

0.04 

thin 

HD117105 

5889 

4.41 

1.13 

-0.29 

3.0 

0.7 

-0.12 

thin 


Notes. The full table is available at the CDS. 


Table 3. Oxygen abundances from the [OI] 6300A line 


Star 

Te// 

(K) 

log(g) 
(cm s' 2 ) 

ft 

(km s' 2 ) 

[Fe/H] 

[Ni/H]" 

EW, Vl 

(mA) 

EW [0 /] 

(mA) 

<5EW[o/] 

(mA) 

[O/H] [oi] 

Gal.Pop. 

Sun 

5777 

4.44 

1.00 

0.00 

0.00 

2.0 

3.4 

0.3 

0.00 

thin 

HD48611 

5337 

4.51 

0.69 

-0.36 

-0.37 

1.2 

1.9 

0.8 

-0.42 

thin 

HD70889 

6051 

4.49 

1.13 

0.11 

0.09 

1.8 

3.3 

0.4 

0.10 

thin 

HD21161 

5923 

4.24 

1.14 

0.09 

0.11 

2.2 

5.9 

2.9 

0.23 

thin 

HD 130322 

5365 

4.37 

0.90 

-0.02 

-0.03 

2.5 

5.1 

1.0 

0.09 

thin 

HD222422 

5475 

4.46 

0.73 

-0.12 

-0.15 

1.8 

2.8 

0.8 

-0.16 

thin 


Notes. (n) [Ni/H] from lAdibekvan et ahl d2012l ). The full table is available at the CDS. 


T e ff at temperatures less than 5600K, for which lower [O/H] 
abundances are found. Although restricted to the low tempera¬ 
ture end, a trend with opposite sign can also be noticeable for 
[06i5s/H]. Nevertheless, if we only consider stars with temper¬ 
atures larger than -5400 K (e.g. when the 6158A line is strong 
enough for reliable measurements) we will find no temperature 
dependence for the 6158A line. In general, given the high abun¬ 
dance scatter observed in these graphs, the small dependence on 
T e ff for cool stars can be neglected. No gravity effect was found. 


3.1. Uncertainties 

There are several sources of random error involved in the pro¬ 
cess of measuring chemical abundances. These are related either 
to the measurement of the EW or to the calculation of the abun¬ 
dances. It is not easy to define the contribution of each error and 
account analytically for them. In the present work we have con¬ 
sidered three different contributions to the uncertainties in our 
measurements. 

First, we studied the sensitivity of the abundances to the stel¬ 
lar parameters: effective temperature T e ff, metallicity [Fe/H], 
surface gravity log(g) and microturbulence //. Variations of these 
parameters introduce small changes in the abundances which, 
regardless of the parameter we vary, also depend on the T e ff 
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Table 4. Average stellar parameters for the three sub-samples of differ¬ 
ent T e ff. 



T eff 

(K) 

log(g) 

[Fe/H] 

ft 

(km s r ) 

Low T ef f 

5292 

4.35 

-0.16 

0.73 

Solar T e ff 

5774 

4.40 

-0.12 

1.00 

High T eff 

6340 

4.50 

-0.10 

1.59 


of the star dAdibekvan et al .112012 ). Therefore, we will carry on 
our analysis with three different sub-samples based on the T e ff 
T < 5377, 5377 < T < 6177 and T > 6177. The middle group 
comprises the stars with T e jf =T 0 ± 400 K. We expect that 
the sensitivity to the stellar parameters behave similar within 
each range of temperature. The average values of the stellar pa¬ 
rameters for each group are shown in Table [4] These parame¬ 
ters together with the average oxygen abundance on each group 
are used to compute theoretical EWs of oxygen lines using the 
ewfind driver of MOOG. Keeping the EW constant, we then cal¬ 
culate abundances by varying, one by one, the stellar parameters 
by an amount equal to their uncertainty (Table 0. Among all the 
values, we select as the final error the largest deviation from the 
original oxygen value. 6158A appears to be more sensitive to 
T e ff while 6300A larger variations with log(g). 

Another source of uncertainty is the statistical photometric 
error due to the noise in each pixel. We can evaluate the partial 
derivative of the EW with respect to the flux dBohlin et al.l ll983: 
ICavrell 19881) to obtain an analytical solution to this contribution: 


EW = 


,6EW = 



2.45e VcrAT 


°< 

E 





Continuum 
Photometric 
Stellar Parameters 


1500 


1800 


2100 


al S/N 


Fig. 6. Contribution to the total equivalent width error of 016158A of 
each source of uncertainty. 


and the ±e continuum placements, where e is the inverse of the 
S/N per resolution element. The width at the continuum level is 
again equal to 6cr. All together yields an error of ±6<xe. 

Figure [6] shows significance of different contributions to the 
total uncertainty of the EW, which is similar for both oxygen 
lines. We have converted the error in abundance due to uncertain¬ 
ties in stellar parameters to the corresponding EW error using the 
ewfind driver of MOOG synthesis code. Each contribution to the 
EW is finally added quadratically to get the one-cr uncertainty. 
Because of the non-linear dependence of the abundance on the 
EW, asymmetric errors will arise when the EW error is propa¬ 
gated to get abundance uncertainties. For further details on the 
determination of the abundance uncertainty we refer the reader 
to the appendix. 


Where C; is the continuum flux level, /, is the flux of each 
measurement, AT represents the wavelength increment per pixel 
( 0.01A for HARPS spectra), e is the r.m.s. relative photomet¬ 
ric accuracy in the continuum and cr is the standard width over 
which we sum six times assuming a gaussian profile. We deter¬ 
mine the standard width theoretically as the quadratic sum of 
different broadenings: instrumental (wavelength divided by the 
spectral resolution of the instrument), macroturbul ence based on 
the relation provided by IValenti & Fischerl ( 2005 :1. microturbu¬ 
lence, rotation (using an average value of 3 kms_') and thermal 
broadening. The latter is calculated with (Gravill992[) : 

A l2kT n 

A A d = — \l -= 4.301 • 10^ 7 T 

c V m 

Where k is Boltzmann’s constant, T represents the effective tem¬ 
perature and p is the atomic weight in atomic mass units. 

The last and the greatest contribution to the uncertainties is 
related to the measurement of the EW and its sensitivity to the 
continuum placement. Although this error can reach a large per¬ 
centage of the EW, it is commonly ignored by many authors, re¬ 
sulting in significantly underestimated uncertainties in the abun¬ 
dances. In order to quantify this contribution, one should cal¬ 
culate the difference between the areas of the gaussian fits by 
placing the continuum in different positions. Given the high S/N 
of our spectra, an therefore the low uncertainty in the rms, the 
possible continuum placements will yield small variations be¬ 
tween the respective Gaussian depths. Thus, the error in EW can 
be approximated as the area of the rectangle between the best 


4. Comparison between indicators 

Abundances from both the OI6158A and [OI]6300A lines were 
measured in 447 stars from the initial sample of 762. This large 
number of stars covers a wide range of temperatures and metal- 
licities, which allows us to carry out a comprehensive compari¬ 
son of the performance of these spectral lines as oxygen abun¬ 
dance indicators. In this section we study the agreement between 
chemical abundances obtained from both spectral lines, and we 
also investigate their reliability and usability related to the stellar 
parameters. 

In section [3 we have already noted that the absolute oxy¬ 
gen abundances from the two lines in the Sun disagree, since 
log e(06i58) is 0.06 dex higher than log HO^oo). The com¬ 
parison for all the stars with common measurements is shown 
in Figure [7] The OI6158A line provides higher abundances in 
stars with [O/H] similar or lower than the Sun, while [OI]6300A 
yields higher values for oxygen-rich stars. This small trend could 
be due to overestimation of the equivalent width when the lines 
are small. However, the trend is not significant as it falls within 
the standard deviation of ±lcr = 0.127 dex. On average, abso¬ 
lute oxygen abundances derived from the 6158A spectral line 
are 0.024 dex higher than [OI]6300A. This value is considerably 
smaller than the typical uncertainties of the oxygen abundances. 
No conclusion can be drawn, except that on average both spec¬ 
tral lines yield compatible values of oxygen abundance within 
the uncertainties. 

As we move towards stars with better quality spectra, i.e. 
higher S/N, the abundance uncertainties decrease. As expected. 
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Table 5. Abundance sensitivities to the stellar parameters of both oxygen indicators for the three sub-samples of different T e /f. 




AT 

eff 

Alog(g) 

A[Fe/H] 


< 



+30 K 

-30 K 

+0.06 dex 

-0.06 dex 

+0.03 dex 

-0.03 dex 

+0.08 kms 

1 -0.08 kms 1 

Low T eff 

A 016158 

-0.03 

+0.04 

+0.02 

-0.01 

0.00 

+0.01 

0.00 

+0.01 

A [OI]6300 

+0.01 

0.00 

+0.03 

-0.02 

+0.02 

-0.01 

0.00 

0.00 

Solar T e ff 

A 016158 

-0.02 

+0.03 

+0.02 

-0.01 

0.00 

+0.01 

0.00 

0.00 

A [OI]6300 

+0.01 

0.00 

+0.03 

-0.02 

+0.01 

-0.01 

0.00 

0.00 

High T eff 

A 016158 

-0.02 

+0.02 

+0.02 

-0.02 

0.00 

+0.01 

0.00 

0.00 

A [OI]6300 

+0.01 

-0.01 

+0.02 

-0.02 

+0.01 

-0.01 

0.00 

0.00 
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Fig. 7. Comparison of the absolute oxygen abundances from OI6158A 
with those obtained from [OI]6300 A lines. Solid lines represent the 
diagonal and the linear fit to the data. Dashed lines enclose 68% of the 
data around the mean (0.024 dex). 



Fig. 8. Difference between absolute oxygen abundances derived from 
the OI6158A and [OIJ6300A lines against nominal signal-to-noise ratio. 
Weighted moving average (boxcar average) of data points is shown as 
a solid line, together with the 1, 2, and 3 sigma dispersion regions in 
different grey shades. Dashed lines represent the simple moving average 
of the abundance errors for each line. In the top panel the Sun has been 
shifted towards lower S/N for plotting purposes. The lower panel shows 
the width of the ± lcr region for the data average and the uncertainties. 


our measurements become more precise and we achieve a 
smaller dispersion in the difference between oxygen indicators. 
Specifically, while the standard deviation for the whole sam¬ 
ple was 0.127 dex, at high S/N it drops to 0.072 dex. How¬ 


ever, for stars with S/N>1000 the average difference in oxygen 
abundances is 0.051 dex, which is compatible with previously 
measured differences in the Sun dCaffau et alJ12008 ). These de¬ 
viations increase as we consider higher S/N ratios, indicating 
a possible discordance between both indicators. The agreement 
between abundance indicators as a function of the S/N is pre¬ 
sented in figure [8] We have performed a weighted moving av¬ 
erage (500 steps, step size 50 and weigthed with 1/cr 2 ) that has 
been smoothed by a factor of 0.4. The 1,2 and 3 standard de¬ 
viations levels from the mean, which are represented by shaded 
areas in the figure, clearly show how the dispersion decreases 
at better S/N, while the mean is shifted towards higher values of 
06158A. This result may suggest that the 0.06 dex disagreement 
found in the Sun is not unique, and may have some physical ex¬ 
planation. 

The different sensitivities of the two oxygen lines to stel¬ 
lar parameters allows us to define windows where each of these 
lines can be used more or less easily. The equivalent width of 
both spectral lines is smaller than 8mA at T e ff >5800K. How¬ 
ever, because of the high excitation energy of 6158A, this line 
becomes significantly stronger in hotter stars. Specifically, the 
EW of OI6158A is always larger than the EW of [OI]6300A in 
all the stars with T e ff >6200K. 



[Fe/H] 


Fig. 9. Fraction of oxygen in the equivalent width of the 6300 A feature 
as a function of metallicity. Dependence on the effective temperature is 
shown in colours. 


It is also important to note the uncertainties that affect each 
line. Because of the density of spectral features in the range of 
6158A, the largest source of error in the equivalent width mea¬ 
surement is the continuum placement. Therefore, it is more com¬ 
plicated to measure this line in cool stars, where in addition to 
many blends, the line itself is very weak. The 6300A feature 
can be considered as problematic in solar-type stars for two im- 


Article number, page 7 ofnu 
















A&A proofs: manuscript no. oxygen 


portant reasons. The first is related to the presence of numer¬ 
ous telluric lines often severely blended with the oxygen fea¬ 
ture. Even more serio us is the problem with the Ni absorption 
line (see above). Ijohansson et al.l ( 2003 ) provided a logfg/') value 
with 14% uncertainty. Our calculations show that the Ni line can 
account for up to 60% of the EW of the 6300A absorption fea¬ 
ture in cool, metal-rich stars (see Fig |9}. Thus, the abundance 
of Ni and its atomic parameters plays a key role in this topic. 
As an illustration, the error in the Ni log (gf) is introducing an 
uncertainty in the oxygen EW measurement of HD 142 equal to 
3.1*qi mA, which is translated into [O/H]=0.2*qqq dex. In stars 
with medium to low S/N this contribution will be smaller than 
the two main sources of uncertainty considered in this work. 
Even though, the larger scatter found in figure [4] with 6300A 
could be partially explained by the uncertainty in the Ni logfg/'), 
in combination with the possible presence of undetected telluric 
lines slightly above the level of the noise. 

The main and the most important aspect of the 6158A line 
is that it does not suffer from known blends. It is as strong as 
the forbidden line in solar analogues and stronger than 6300A 
at temperatures higher than 6200K. Our results suggest that the 
6158A line can be used as a reliable indicator of oxygen abun¬ 
dance in solar-type stars. However, we should bear in mind 
tha t it may s uffer non-negligible non-LTE effects in hot stars 
(iSitnova et al.ll2013b . 

5. Oxygen trends in the thin and thick disks 

5.1. Observed dispersion of [O/H] 

The abundance scatter seen in the Galactic trends of many chem¬ 
ical elements is very small. [Adibekvan et al.l d2012l) carried out 
a homogeneous study and found small abundance dispersion in 
such elements as Si or Ti. The abundances presented here have 
been derived from the same spectra and stellar parameters. Thus, 
using the same database and techniques as [Adibekvan et alJ 
( 2012 ) we would expect to get equally precise and trustworthy 
results for other elements. However, those elements have many 
more reliable spectral lines for abundance analysis, and there¬ 
fore is not surprising that we found a wider scatter in the case 
of oxygen (fig. HI. Even so, scatter found in the present work for 
[C> 6 i 58 /H] is compar able to that of Mg, whic h was derived using 
three spectral lines (Adibek van et al] [20121) . IPetigura & M arcv 
(1201 1 1) suggested an astrophysical origin to explain the disper¬ 
sion found in oxygen trends. 

We consider a sub-sample of 296 stars, 55% of the sample, 
for which we have errors smaller than 25% in EW in order to ver¬ 
ify whether the observed abundance scatter of oxygen is due to 
the lack of precise measurements. We also restricted the analysis 
to the OI 6158 line, since we found in this study that it is more 
trustworthy. We see in Fig. |T0]that the scatter is still very high 
in this ‘most precise’ sub-sample at [Fe/H] < -0.2 dex. On the 
other hand, at higher metallicities the scatter is not remarkably 
reduced by selecting the best measurements, and only outliers 
are removed. Error bars are about the size of the dispersion and 
it is therefore difficult to draw any conclusion about the origin of 
the scatter. 

The observed scatter at a given metallicity is expected to 
have two contributions: a real star-to-star scatter, and also a dis¬ 
persion due to the uncertainties in our measurements. We have 
performed a moving average (250 steps with 0.03 dex width) to 
derive both the average [O/H] with their respective scatter, and 
the mean +<5[0/H] as a function of metallicity. The goal of this 



Fig. 10. [O/H] vs. [Fe/H] for OI 6158 A. Filled (open) circles represent 
data for which EW measurements have an uncertainty under (above) 
25%. Absolute error bars are over-plotted for our best measurements. 
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[Fe/H] 

Fig. 11. Upper panel: [O/H] vs. [Fe/H] for OI 6158 A. Only stars with 
EW relative error below 25% are considered. Solid line represent the 
average oxygen abundance. The 1,2 and 3 sigma levels of dispersion 
around the mean are presented as shaded areas. Dashed line shows the 
average uncertainties. For plotting purposes the average uncertainties 
are shown around the mean abundance. Lower panel: Difference be¬ 
tween the ±1 <x width of the scatter around the mean and the ± l<x width 
of the uncertainties. 


test is to find out which of the two contributions to the scatter 
dominates our sample. Results of the test are depicted in fui.lTil 
We are aware that this test is not valid at the edges of the metal¬ 
licity range due to the undersampling of the bins, which makes 
unreliable the provided mean values. 

In order to evaluate the relative importance of the contribu¬ 
tions to the scatter, in the lower panel of the figure we show 
the difference between the widths of the ±l<x regions. A nega¬ 
tive residual would imply an overestimation of our uncertainties, 
while a positive residual could be due to either a real scatter or 
to an underestimation of the errors. Since the uncertainties have 
been derived with the same method for all stars, and there is 
no dependency on metallicity, differences in residuals as a func¬ 
tion of metallicity are expected to be caused by real star-to-star 
dispersion. At high metallicities the observed scatter and the un¬ 
certainties are comparable, and do not allow us to detect the pos- 
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Fig. 12. [O/Fe] vs. [Fe/H] for [01] 6300 A and 01 6158 A. Stars from different populations are separated according to their chemical composition. 




sible presence of an astrophysical scatter. However, the positive 
residual found at [Fe/H]<-0.2 dex reveals a real star-to-star scat¬ 
ter that cannot be explained by uncertainties. The results of this 
test are not surprising, and they are in perfect agreement with the 
fact that at [Fe/H] <-0.2 dex there is a mixture of stars from two 
different populations of the Galactic disk (see below for details). 
However, this expected results confirm the effectiveness of this 
method to show up the presence of a real scatter. In addition, the 
satisfactory agreement between uncertainties and scatter at high 
metallicities support the validity of our uncertainty calculations. 
Given that the test was performed using a sub-sample of the best 
measurements, this last result should consequently be interpreted 
with caution. A similar test with the whole sample would show 
a slight overestimation of the uncertainties (negative residuals), 
due to the exponential growth of the errors as we move towards 
lower S/N. Therefore, we conclude that although our uncertainty 
calculations are overestimated at very low S/N, they appear to 
satisfactory account for the observed scatter at medium to high 
S/N. 


5.2. Thin and thick disk separation 


It is believed that q elements, such as oxygen, are produced in 
massive stars that enrich the interstellar medium when they ex¬ 
plode as type II supernova. By contrast, iron peak elements are 
mostl y synthesized in SNe la explosions (e.g. lThielemann et al.1 
2002 ). The oxygen over iron ratio can provide hints on the for¬ 
mation of the Galactic disk and its chemical history, but the ob¬ 
served Galactic trends are diffused as a result of the mixture of 
stars from two different populations in the studied metallicity 
regime. 

The chemical distinction between thin and thick disks was 
first noticed by iFuhrmannl (1 19981) and further investigated by 
iProchaska et al.l fl2000 ), among others. We have classified the 
stars in our sample between thin disk, high-q metal-poor (thick 
disk) and hi gh-q metal-rich stars acc ording to chemical criteria 
described in lAdibekvan et ah ( 201 1). These authors show that 
disk populations are clearly separated in terms of [q/Fe], where 
‘a’ refers to the average abundance of Mg, Si, and Ti. The ra¬ 
tio of oxygen over iron for the different disk populations is pre¬ 
sented in Figure[T2] There is no obvious difference seen between 
the thick and the thin populations from [OI]6300, which can be 
attributed to the larger errors and scatter. We see a clear over¬ 
abundance of oxygen in thick disk stars once the OI6158A line 
is considered. Moreover, using the sub-sample with the most pre¬ 



[Fe/H] 


Fig. 13. [O/Fe] vs. [Fe/H] for the stars with 5 EW < 25%. Stars from 
different populations are separated according to their chemical compo¬ 
sition. Solid lines represent the fitted trend, weighted with the uncer¬ 
tainties, for each population (thick and humr stars are fitted as a single 
population). Absolute error bars are over-plotted. 


cise abundance values (Fig. E}. we can check that none of the 
thick disk stars falls below the trend for the thin disk, although 
few thin disk stars have [O/Fe] higher than expected for this pop¬ 
ulation. In addition, both trends appear to be remarkably sepa¬ 
rated at lower metallicities. Thus, we can confirm with our data 
that oxygen, just like the other a elements, is enhanced in the 
thick disk stars. 

The high-q metal-rich pop ulation (hqmr) was presented by 
lAdibekvan etaD (1201 ll [20131) as a new population due to the 
gap in metallicity found between the latter and the thick disk. 
In addition, a gap in the distribution of [q/Fe] at 0.17 was 
al so noticed. The present work is based on the same sample 
as lAdibekvan et alJ (1201 ll ). and therefore our data also show a 
gap in metallicity. The separation in the [O/Fe] distribution is 
not as clear as in [q/Fe] because both populations are mixed 
around [O/Fe]~0.25, probably caused by the larger errors in the 
determination of oxygen abundances compared to other q ele¬ 
ments. Our most precise sample show that all the stars from 
the hqmr population fall above the average trend for thin disk 
(Figfl3l). and thus we can confirm that hqmr stars are enhanced 
in [q/Fel (where q refe rs here to oxygen) . Recent studies by 
iBensbv et al] d2014l) and lHelmi et al] (120141 ) also confirmed the 
alpha enhancement of a metal-rich population, nevertheless, they 
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did not observe a gap in metallicity between them and the thick 
disk stars. iBensbv et al.l (120141) argued that the high-ar metal-rich 
family might be a metal rich tail of the thick disk. With this new 
oxygen abundance determination we cannot draw conclusions 
about this population, i.e. whether is distinct or just a continua¬ 
tion of the thick disk towards higher metallicities. 

Once the disk populations are separated, we can revisit the 
issue about the origin of the scatter in the observed [O/H] trends. 
We repeat the test performed for the best sample, but this time 
independently for the thin and thick (including hrmir) disk pop¬ 
ulations. The moving average for the oxygen abundances and 
+b[0/H] uncertainties are shown in figure [141 together with the 
residuals of the ±lcr width regions. Given the small residuals 
found between uncertainties and scatter for thick and homr pop¬ 
ulation, and taking into account the reduced number of stars that 
we have in this group, no conclusion can be drawn about the 
existence of a real star-to-star dispersion. In contrast, thin disk 
stars seem to show a positive residual at -0.2<[Fe/H]<0.3, which 
reaches 0.03dex at [Fe/H |~0.1. This residual may reveal a real 
astrophysical scatter, since uncertainties are unlikely to be un¬ 
derestimated. On the other hand, it can be checked that the thick 
and thin disk populations appear to be separated at a lcr level, 
thus confirming the enrichment of cr-elements in the thick disk. 
In addition, the enhancement of the hamr population is clearly 
revealed. 
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[Fe/H] 

Fig. 14. Upper panel: [O/H] vs. [Fe/H] for OI 6158 A. Only stars with 
EW relative error below 25% are considered. Stars are sepparated in 
different disk populations according to a chemical criteria. Solid line 
represent the average oxygen abundance. The lcr level of dispersion 
around the mean is presented as a shaded area. Dashed line shows the 
average uncertainties. For plotting purposes the average uncertainties 
are shown around the mean abundance. Lower panel: Difference be¬ 
tween the ± lcr width of the scatter around the mean and the ±1 <t width 
of the uncertainties. 


5.3. Comparision with previous studies 


Several studies have addressed the behaviour of oxygen over 
iron in the thin and thick disk populations. In general, all of 
them agree that thick disk stars show an enhancement of a el¬ 
ements, being this caused by different ratios of SNII/SNIa in 
both disks. Some of them found a shallow decline in the 
trend at high [Fe/Hj for thin disk stars (e.g. Bensby et al 


Rami rez et al J 12007 : iTakeda & Hondal[2005 


Petigura & 


O/Fe] 

20041 


Vlarcvl 


1201 It IBensbv et al.l l2014h . indicating that in the thin disk both 
SN la and SN II contribute at a steady rate to the enrichment 
of the interst e llar medium . In contrast. [Ramfrez et al.l(l2013l) and 
iNissen et al.l (120021120141) results suggest a flattening of [O/Fe] 
at super-solar [Fe/H]. Our data also suggest that flattening in 
the [O/Fe] ratio may appear at solar metallicities (Fig.lT2l>. The 
“knee” is inappreciable in the case of [OI], and thus, for the mo¬ 
ment we cannot rule out that a change of slope below [Fe/H] ~ 
0 is real. This apparent lack of agreement between both oxygen 
indicators at high [Fe/H] can be attributed either to the scatter 
in the [OI]6300 data or because some blend may appear at high 
metallicities near the OI6158A line, leading to an overestimation 
of the EW. However, although further work needs to be done to 
study the linelist near 6158A spectral line, neither the FWHM 
or the symmetry of the line seem to suggest the presence of an 
unknown blend. 


In order to study the behaviour of [O/Fe] towards the metal- 
poor regime in Fig. |T5] we compare our results of OI6158A 
with those of lRamfrez et al .i (12013D . These authors derive oxygen 
abundances from the infrared triplet and correct them for non- 
LTE effects, but neglect inelastic collisions ^with hydrogen. We 
also make a comparison with the results of Bensbv et al.l (l2014l> 
and INissen et al.l (f2014 ). both of them based on a combination 
of samples for which oxygen was measured using [OI]6300 
and the OI IR triplet, and only the OI IR triplet, respectively. 
INissen et al.l (12014 1 corrected oxygen abundances for departures 
from LTE accounting for H collisions, as described by the clas- 
sical Drawin formu la dDrawi nS ll968l h and scaled by S//=0.85. 
Bensbv et al. (2014) used for this purpose an empirical correc¬ 
tion derived by comparing abundances of 60 stars from [OI]6300 
and 01774 indicators. 


Wh ile our results mostly a gree with those of IBensbv et alJ 
l20fi and N issen et al] (120141), we find clea r disagreement for 
the metal-poor stars of lRamfrez et al.l (1 201 3l) . Fig. Q3] shows at 
[Fe/H] <-0.5 how the [O/Fe] rise becomes steeper and we obtain 
much larger values for [O/Fe] than lRamfrez et al] (l2013i) . On the 
other hand, despit e the similar slope , our results reach higher 
[O/Fe] ratios than INissen et alJ (l2014l> and IBensbv et al] {2014} 
at [Fe/H]—0.8 dex. Although the most metal-poor stars in our 
sample have weak OI6158A oxygen lines, with EW between 2 


and 4 mA, we can safely rule out the possibility that our values 
are overestimated because of the poor continuum measurements 
and/or other sources of errors. To illustrate this point, in Fig. [T6] 
we show one of the targets in our sample HD 22879 with one of 
the highest [O/Fe] =0.81 and [Fe/H |=—0.82. The combined spec¬ 
trum of this star has S/N=1000 and our precise measurements 


of OI6158A EW, which perfectly match the synthetic spectra, 
clearly rule out an EW weaker than 2.5 mA. Thus, our best value 
for this star is [O/Fe]=0.81. These high oxygen over iron ratios 
are also found with [OI] measurements, as an example HD25704 
is a star with [Fe/H] =-0.83 and [O/Fe] =0.69 dex. Finally, it is 
noteworthy to mention that the observed trend towards higher 
values of [Of,i.WH I at low temperatures (Fig. 0) do not affect the 
observed [O/Fe] trend since none of the high [O/Fe] stars have 


T eff <5600K. 

Given the quality of the observations and measurements pre¬ 
sented in these studies, the simplest way to explain the dis¬ 
crepancy with other works rely on the different non-LTE cor- 
rection s ap plied to the abundance calculations. Differences with 
iRamfrez et alJ ( 2013 ) arise because the inelastic collisions with 
H atoms were neglected in their OI IR triplet calculations, which 
would have shifted the trend towards the LTE results. Actually, 
calculations of non-LTE bv lFabbian et al.l ( 2009 ) taking into ac- 
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Fig. 15. Comparison of our [CLiss/Te] vs. [Fe/H] trends to those derive d by Ramirez et all(|20l3l). IBensbv et al.l ll20 1 4l) an d lNissen et al] d2014tl . 
For comparison p urposes, abundance v alues from Ramirez et al. (2013) and Nissen et al. (2014) have been rescaled to the same solar reference 
(log f (Q) ,=8.71). IBensbv et al] d2014h performed a normalization relative to their solar values on a line-by-line basis. However, individual line 
abundances are not provided, which makes not possible to rescale to the same solar reference. 


count inelastic collisions with neutral oxygen conclude that non- 
LTE effects presented by these authors were overestimated by 
~0.10 dex at [Fe/H] —1.5 dex. On the other hand, the disagree¬ 
ment found between [OI]6300 and 01IR triplet in the Sun lead 
us to question the non-LTE correction applied by IBensbv et alJ 
( 2014 1). Although this correction may be a good first approach 
to the problem, it would be useful to determine it with more ac¬ 
curacy. Given the current knowledge about departures from the 
LTE, and the weight of hydrogen collision to bring the system 
closer to the LTE, iNissen et alJ (120141) approach appears to be 
the most reasonable. 

The effects of departures from LTE for 6158A transitions 
have been widely studied in the Sun. Unfortunately, there is not 
much literature about non-LTE corrections at lower [Fe/H] for 
this spectral line. As a high-excitation line, it is formed in deep 
layers of the atmospheres, where the source function and line 
opacity effects work in opposite directions. Thus, the combined 
effect results in very small non-LTE effects in the Sun —0.03 dex 
(iCaffau et alJ 12008 1). However, generally speaking, departures 
from LTE are more severe towards higher T e ff and lower log(g) 
and [Fe/H]. Therefore, one would expect larger non-LTE correc¬ 
tions in the low-[Fe/H] tail of our [O/Fe] distribution. Tempera¬ 
ture is not expected to increase the effect of non-LTE in our stars, 
for which effective temperatures are between 5600 to 6200K at 
low metallicitv. lSitnova et al . (120131) studied non-LTE effects for 
oxygen permitted transitions in different stars and for a grid of 
parameters, taking into account hydrogen collisions (S//=l). Ac¬ 
cording to their work, non-LTE corrections are not expected to 
be larger than -0.05 dex in our sample. 

6. Discussion 

In this article we show that the [O/Fe] ratio rises steeply with de¬ 
creasing metallicity, and reaches 0.8 around [Fe/H]=-1. Our data 
are not extended below [Fe/H] = -1 and it is not clear if this ratio 
will stay constant or keep rising in halo stars down to [Fe/H]<-3. 
Obviously both oxygen lines employed in this work will disap¬ 
pear in halo dwarfs since their EWs will be less than 1mA. Thus, 
there is no chance to extend this study to F type metal-poor halo 
dwarfs. Finally, we note that despite considerable observational 
effort the trend of the [O/Fel ratio in the halo is still uncertain . 
However, most of the studies (llsraelian et al.l200lllTakedaf2 003: 
Fulbright & Johnson 2003) suggest that main-sequence dwarfs 
provide more reliable and consistent abundances than giants. 

A steep rise of the [O/Fe] ratio and high values above 
0.7 found in the present analysis can hardly be explained in 
the framework of classical models of galactic chemical evo¬ 


lution. We consider, for example, the models proposed by 
IChiappini et al.ld2001l) . This model considers two main accretion 
episodes for the formation of the Galaxy, the first forming the 
halo and bulge in a short timescale and the second one forming 
the thin disk. They propose a short halo formation timescale of 
about 0.8 Gyr. At a certain epoch (coinciding with the halo-disk 
transition) core-collapse SNe, responsible for the production of 
O, stopped exploding while Fe, produced by the long-living SNe 
la, continued to be produced. Star formation must have ceased 
for a period that cannot be longer than 1 Gyr. The two-infall 
model fits the observed metallicity distribution of the G-dwarfs 
by assuming a long timescale for the thin-disk formation. These 
models have been calculated assuming a constant plateau at 
[Fe/H]<-1 for all a-elements including oxygen, for which they 
propose [O/Fe]=0.5 at [Fe/H ]=-1. These models predict a light 
increase in [O/Fe] towards lower [Fe/H] reaching rO/Fe]~0.4 a t 
[Fe/H]=-1. However, some authors (e.g. I Boesgaard et alj|201 In 
have already reported that a-elements do not show a plateau 
at [Fe/H]<-1. Moreover, the steep rise of [O/Fe] does not nec- 
essarily imply the same behaviour for other a-elements (e.g. 
[Ramatv et al.l 120001) . The different trend for oxygen compared 
to the other a-elements at [Fe/ H]>0 can also be explained by 
the fact that oxygen is only produced in Type II SN. The lev¬ 
elling out of the [a/Fe] trend for other a elements is clearly ex¬ 
pected since they also have small contributions from Type la SN. 
Metallicity dependent oxygen yields will make the overall pic¬ 
ture more complete and may provide additional clues to explain 
these observations dPrantzos et all 19941) . 

Any model for the thick disk origin should account for sev¬ 
eral observational facts. The first one is the old age of thick disk 
stars and the temporal gap (1-2 Gyr) with respect to the old¬ 
est thin disk stars. Then, the tight relationship found for sev¬ 
eral abundance ratios in the thick disk. And finally, the varia¬ 
tion of the abundance ratios with metallicity. It will be an in¬ 
teresting challenge for these models to explain observations of 
a-elements from the same HARPS database dNeves et al.ll2Q09i: 
Adibekvanet ail l2012l) as well as [O/Fe] trend ratio presented 
here since they come from the same homogenous sample. 

Different models have been proposed to explain observa¬ 
tional characteristics of the thick disk. It is possible that thick 
disks have been c r eated through accr etion of galaxy satellites 
dMeza et al.l 120051 : lAbadi et~akf l2003h . where thick-disk stars 
therefore have an extragalactic origin. It has also been pro¬ 
posed that thick disks were born thick at high redshift from 
the internal gravitational instabilities in ga s-rich clumpy disks 
(iBournaud et al] l2009t iForbes et al] 120121) . Another possibil¬ 
ity is that thick disks are created through the heating of pre- 
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Fig. 16. Spectral synthesis of the oxygen line at 6158A in HD 22879 
(upper panel) and [OI]6300 in HD 25704 (lower panel). The observed 
spectra (dots) with the lcr photometric noise (rms) is compared with 
the synthetic spectra, computed for the corresponding [Fe/H] and 
[O/H] abundances obtained from our EW measurements (solid line). 
Dashed lines represent the synthetic spectra for the ±1 <x error abun¬ 
dances, which are (5[Fe/H]= ±0.03 dex, 8[0^i/H] dex and 

S\0 63K) /H\ =^q Jg dex (propagated from the provided EW error). Also 
shown are the stellar parameters, nominal S/N, abundances. EW and 
residual intensities (observed - predicted). 

existing thin disks with th e help of mergers (lOuinn et alJI 19931: 
I Villalobos & Helmilt2008l) . whose rate decreases with decreas¬ 
ing redshift. One could propose that a massive merger will cre¬ 
ate very intensive star formation at [Fe/H]=-1 and produce stars 
with very high [O/Fe] ratios. 

The thick and thin disks have formed on different timescales. 
Our data shows that the average ages of thin and thick disks are 
4.73 Gyr and 9.20 Gyr, respectively. Details about age determi¬ 
nation are provided by iDelgado Mena et al.l (!2014t) . In fact, all 
our stars with [O/Fe] >0.6 belong to the thick disk and are older 
than 6 Gyr. The Galactic thin disk has not had such an intense 
star formation history as the thick disk. The shallow decline in 
the [O/Fe] is caused by continuous star formation with no fast 
initial enrichment from Type II SN. 

In the recent past there has been a growing conviction that 
radial migration can be responsible for the formation of thick 
disks by bringing out high-velocity-dispersion stars from the in¬ 


ner disk and the bulge (e.g. Sch onrich & Binnevll2009l) . While 
by now it has been established that because of conservation of 
vertical action migra t ion does not contribute to disk thickenin g 

dMinchev et al.ll2012t iMartig et alJl2014l IVera-Ciro et al.ll2014l) . 

radial migration could bring stars with high [O/Fe] from inner 
galaxy to solar neighbourhood dMinchev et alj|2013l) . 

The stellar ev o lution calculations for massi ve stars 
(IWooslev & Weaverl 1 1 995b iThielemann et al] 1 19961) predict 
that the ejected mass of (O/Fe) by a massive star is an increasing 
function of the stellar mass. Thus, an increase of the slope of the 
IMF during the halo/thick disk transition phase could possibly 
explain our observations. A very steep linear increase of the 
[O/Fe] in thick disk stars could also be due either to a lower 
contribution of Fe from Type II SNe or to a metallicity depen¬ 
dent yield of oxygen (or a combination of these two effects). 
The change of metallicity on oxygen yield has a different impact 
depending on the mass range of stars considered (Maederi 1992]) . 

Perhaps the easiest way to explain high [O/Fe] ratios in 
metal-poor stars is by modifying the so called mass-cut parame¬ 
ter in the models of explosive supernova. Present explosive nu¬ 
cleosynthesis calculations for Type II SN are based on induced 
supernova explosions by either depositing thermal energy or in¬ 
voking a piston with a given kinetic energy. Induced calcula¬ 
tions (lacking self consistency) utilize the constraint of requiring 
ejected 56 Ni masses from the innermost explosive Si burning lay¬ 
ers in agreement with supernova light curves being powered by 
the decay chain 56 Ni- 56 Co- 56 Fe. This can also serve as guidance 
to the supernova mechanism with mass cuts based on 56 Ni in the 
ejecta. It means that the position of a mass cut can define the 
amount of Fe eject ed from the supernova as well as the [O/Fel 
ratio in the ejecta (Wooslev & Weaver 1995; Thielemann et al.l 
1 19961) . 

7. Summary and conclusions 

We present a uniform study of oxygen abundances in a large 
sample of F- and G-type dwarfs from the HARPS database. A 
homogeneous set of the atmospheric parameters was adopted, 
and two independent analyses for different indicators were per¬ 
formed. We provide the first accurate and homogeneous compar¬ 
ison of the OI6158A and [OI]6300A lines of oxygen. 

Oxygen is one of the most controversial elements owing to 
the problems related to inconsistent abundances derived from 
different spectral lines. In this work we have investigated the be¬ 
haviour of the lines at 6158A and 6300A. Our measurements 
suggest that abundances derived from both oxygen lines are 
consistent within their respective error bars. However, we have 
found that for spectra with a very high S/N these values are no 
longer compatible within errors. Mainly, the abundance derived 
from the OI6158 A line is 0.072 dex higher than that derived from 
the [OI]6300A line. A similar difference has been previously re¬ 
ported for the Sun dCaffau et al.ll2008h . Thus, our results indicate 
that the disagreement is not unique, and that there must be some 
physical explanation for this fact. 

For F-type stars with temperatures higher than 6200 K, the 
use of OI6158A is recommended. The equivalent widths of this 
line are larger than those of the [OI] 6300. Thus, the errors due 
to the placement of continuum are lower for OI 6158. This line 
can also be useful in the super-solar metallicity regime where the 
Ni 6300 line may contribute up to 50% of the feature at 6300A. 
In addition, the uncertainty of log( gf) for Ni is non-negligible; 
therefore, errors in the measurement of [OI]6300A could be sig¬ 
nificant. In contrast, at low temperatures and metallicity below 
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solar values, OI6158A becomes weak and yields more uncer¬ 
tain values of oxygen abundances. In any case, both lines can 
be used in solar analogues if S/N > 300, allowing us to measure 
equivalent widths as small as 2 mA. With this quality of data, 
the OI6158A line also represents a good alternative when the 
[OI]6300A line is severely blended with telluric features. It is 
always good to measure both lines whenever possible in order to 
check their consistency. 

We found that the abundances derived from the OI6158A 
line suffer from much smaller dispersion when considering the 
trends [O/H] and [O/Fe] vs. [Fe/H], in contrast with those de¬ 
rived from the [OI]6300A line. The scatter found in the trends for 
a sub-sample of 252 high-precision measurements of OI6158A 
is comparable to that o f Mg, derived using three spectral lines 
(Adibekvan et al. 2012). In addition, our results suggest that the 
[O/Fe] ratio reaches a constant value at super-solar metallicities, 
in contrast with some previous results that reported _a shallow 
declin e. A similar trend is see n in the results of Ramire z et al.l 
d2013h and iNissen et al . (20H). 

Although there is no clear boundary between oxygen abun¬ 
dances in the thin and thick disks, the stars studied in this article 
show that some bifurcation might exist at [Fe/H] < -0.3 dex, 
as found with the high-precision sub-sample of OI6158A. Our 
results also indicate a clear enhancement of oxygen abundances 
in the thick disk stars, as has been widely reported for other a- 
elements. 
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Appendix A: Statistical approach to the 
determination of chemical abundances 


It is well known that the relation between the equivalent width 
and the abundance of a spectral line is non-linear. If we assume 
that the EW has a normal probability distribution, a non-linear 
transformation will yield an asymmetric abundance distribution. 
The most obvious consequence of this fact is the asymmetric 
abundance errors that one can find in any abundance analysis. 
However, these asymmetries are negligible for strong spectral 
lines. The present analysis is based on weak lines of oxygen, 
whose EW is smaller than 5 mA in more than 50% of the stars. 
Moreover, EW confidence intervals are not preserved under this 
non-linear transformation, which means that abundance errors 
derived from an LTE analysis does not correspond to the same 
tolerance region. Aiming to provide statistically meaningful un¬ 
certainties, we carry out a detailed statistical error analysis. 

If we assume a normal probability distribution for the equiv¬ 
alent width, truncated at_zero, we can obtain the abundance dis¬ 
tribution as (lBarlowll2003l) : 

G t (EW) 

P([X/H]) = |, ry I, B | 

\d[X/H]/dEW\ 

where Gr(EW) is the adopted distribution for the EW with mean 
equal to the measured EW and standard deviation calculated as 
the final lcr uncertainty of the EW. A simple relation between 
abundance and EW is provided by the curve of growth. Consid¬ 
ering a single spectrum line, we can write: 

mlog(EW) = [X/H] + c 


where m and c are constant coefficients. For each star we can find 
these coefficients by evaluating two values of EW, and deriving 
their abundances following the standard LTE analysis. The prob¬ 
ability density function of the abundance will be: 


P([X/H]) = 


MO 

2<X V27T 


10 - 


.(io“S“) 


We can test the correctness of this model by randomly generating 
a set of EW values following a normal distribution, truncated at 
zero. The assumed abundance distribution should be able to re¬ 
produce the distribution of the random EWs, transformed into 
abundance using the curve of growth. This test is shown in Fig¬ 
ure [AT] 

Because of the non-linearity of the transformation, this dis¬ 
tribution is remarkably asymmetric. Several problems will arise 
from this fact, such as the asymmetry of the abundance errors or 
the biased result towards the high tail which is obtained from the 
LTE abundance analysis. 

Ideally, the normal distribution of the EW is not truncated, 
and the measured EW of the spectral line corresponds to the 
mean, median and mode, as expected for normal distributions. 
Under monotonic transformations the median is preserved and 
therefore the abundance from the LTE analysis will be the me¬ 
dian of the distribution, acceptably quoted as a result, but no 
longer the mean or mode. However, when we address the mea¬ 
surement of small lines, the normal distribution is truncated at 
zero as a result of the negative values of EW. In this situation 
the measured EW is only the mode of the distribution (not the 
median), and thus its transformation will not yield the median 
of the abundance probability distribution: the abundance derived 
in the LTE analysis will not only be biased towards the left tail 
of the distribution, but will also be smaller than the median. In 



[X/H] 


Fig. A.l. Abundance distribution probability for a large asymmetry 
case. Red dots represent the analytically derived function. Green bars 
represent the distribution of the abundance calculated with the curve of 
growth for a set of random EW values that follow a truncated gaussian 
distribution. 


this case, it becomes especially useful to calculate the probabil¬ 
ity distribution in order to provide statistically consistent results. 
Abundances derived from the LTE analysis have been used in 
the present work; however, we provide in Tables IA.2I and IA. II 
the statistical results. Figure [A~2l shows the deviation of the LTE 
abundances from the median value of the statistical distribution. 
It is shown that for EWs larger than 8 mA , regardless of the 
signal-to-noise, the result from the LTE analysis coincide with 
the median of the abundance statistical probability distribution. 

Abundance uncertainties are often calculated separately and 
added together in quadrature (positive and negative contribu¬ 
tions independently). Nevertheless, this procedure is not justi¬ 
fied mathematically (lBarlowll2003l) . Throughout this paper we 
provided the errors in EW, which are symmetric. However, the 
uncertainties in abundance should be calculated using the prob¬ 
ability distribution. The errors are calculated as the abundance 
values where the probabilities are 0.159 and 0.841 for the nega¬ 
tive and positive errors, respectively. These abundance intervals 
would enclose the +1 <x region around the median, in a normal 
distribution. We present in Fig. lA.3l the median abundance with 
the corresponding uncertainties. Tables [A2l and I aTI show abun¬ 
dance results from this statistical analysis for the OI6158A and 
[OI]6300A lines, respectively. 
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Fig. A.2. Deviation of the LTE abundances from the median value of the statistical probability distribution. Dependence on the nominal signal-to- 
noise is shown in colours. 




Fig. A.3. [O/H] vs. [Fe/H] for both oxygen abundance indicators, where [O/H] is taken as the median of the statistical distribution. The yellow dot 
corresponds to the solar abundances. [O/H] absolute error bars are over-plotted. 


Table A.l. Statistical oxygen abundances from the OI 6158 A line 


Star 

EW 0/ (mA) 

(5EW 0 / (mA) 

[O/H ] 0 i 

Median [0/H] 0/ 

cr oi 

+<S[0/H] 0/ 

-<S[0/H] 0/ 

Sun 

3.7 

0.3 

0.00 

0.00 

0.04 

0.04 

-0.05 

HD70889 

6.9 

0.5 

0.10 

0.10 

0.04 

0.04 

-0.04 

HD21161 

7.5 

2.8 

0.18 

0.18 

0.27 

0.18 

-0.26 

HD37226 

7.4 

0.8 

-0.01 

-0.01 

0.06 

0.06 

-0.06 


Notes. The full table is available at the CDS. 


Table A.2. Statistical oxygen abundances from the [OI] 6300 A line 


Star 

EW [oi] (mA) 

<5EW[oj] (mA) 

[O/H] [o/] 

Median [0/H] [0 /] 

C[0/] 

+(5[0/H] [0 /] 

-<5[0/H] [0/J 

Sun 

3.4 

0.3 

0.00 

0.00 

0.04 

0.04 

-0.05 

HD48611 

1.9 

0.8 

-0.42 

-0.42 

0.26 

0.16 

-0.25 

HD70889 

3.3 

0.4 

0.10 

0.10 

0.05 

0.05 

-0.06 

HD21161 

5.9 

2.9 

0.23 

0.24 

0.29 

0.17 

-0.28 


Notes. The full table is available at the CDS. 
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